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We study the spinless and spinful extended Hubbard models with repulsive interactions on the 
kagome and the decorated honeycomb ("star") lattice. Using Hartree-Fock mean-field theory, we 
show that interaction-driven insulating phases with non-trivial topological invariants (Chern num- 
ber or Z2 invariant) exist for an experimentally reasonable range of parameters. These phases 
occur at filling fractions which involve either Dirac points or quadratic band crossing points in 
the non-interacting limit. We present comprehensive mean-field phase diagrams for these lattices 
and discuss the competition between topologically non-trivial phases and numerous other ordered 
states, including various charge, spin, and bond orderings. Our results suggest that Z2 topologi- 
cal insulators should be found in a number of systems with either little or no intrinsic spin-orbit 
coupling. 
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I. INTRODUCTION 



The study of topological properties of a quantum sys- 
tem with many degrees of freedom can provide insights 
into global features of ground states and can reveal phys- 
ical behaviors which are robust against weak perturba- 
tions. While the terminology of topological order has been 
used to describe different aspects of a quantum system 
with interactions,^^'' we will focus on topological prop- 
erties which are stored in the set of single-particle wave 
functions describing band structures of materials with a 
bulk gap. One famous example is the integer quantum 
Hall effect where the topological property is encapsulated 
in an integer called the Chern number. It has been shown 
that the Chern number is directly related to quantized 
values of the Hall conductivity."* A nonzero Chern num- 
ber requires the breaking of time-reversal symmetry ei- 
ther by an external magnetic field or, in the absent of a 
net magnetic flux through a unit cell, by microscopically 
circling currents.'^ 

Breaking of time-reversal symmetry is not necessar- 
ily required to define topological invariants which dis- 
tinguish different bulk insulators (or superconductors). 
Based on the random matrix theory, a comprehen- 
sive classification scheme for non-interacting systems has 
been worked out.^ Among all classifications, topologi- 
cal insulators (TIs) with time-reversal symmetry have 
raised considerable interest in recent years (see Refs. [7- 
9]). TIs are well described by conventional band the- 
ory. However, they are a distinct phase of matter with 
bulk energy gaps and an odd number of time-reversal 
symmetry protected gapless modes on their edge (sur- 
face in three dimensions). ^°^*^ In two dimensions, it is 
also termed the quantum spin Hall state (QSH). This 
state is distinct from ordinary insulators by a nonzero 
value of a Z2 invariant. ^"'^* In three dimensions, there 
are four Z2 invariants characterizing either a strong topo- 



logical insulator, a weak topological insulator, or a triv- 
ial insulator. ^^"*'' The Z2 invariants can be obtained via 
knowledge of the single-particle wave functions alone. 

The key to experimental realizations of TIs (at least 
so far) is strong intrinsic spin-orbit interaction originat- 
ing from relativistic effects. The topologically nontrivial 
behavior in these systems is stabilized by a strong spin- 
orbit coupling which leads to a "band inversion". ^^"^^ 
While the experimental search for the TIs in real ma- 
terials with strong spin-orbit coupling is still under way 
with a number of examples found to date,**^^^ the cur- 
rent theoretical research in TIs is quite diverse. On the 
one hand, there have been intensive first-principle stud- 
ies to identify potential candidate materials for TIs.^^^'^" 
On the other hand, the study of TIs in the presence 
of disorder''*'''^ and interplay of spin-orbit coupling and 
electron-electron interaction'^'''''^ have been carried out. 
New exotic phases have been proposed, such as a topolog- 
ical Mott insulator,'^^ which has a gapped charge sector 
but gapless spinon excitations on the boundary. 

In the present paper, we focus on yet another class 
of systems in which the topologically nontrivial na- 
ture of the wave functions is a result of spontaneously 
broken symmetry in an interacting system. ■''"'"'''' These 
interaction- driven topological insulators possess conven- 
tional order parameters and the topological order is 
locked to those. Microscopically, the topological phases 
are described by the spontaneous generation of (spin) 
currents, a popular theoretical idea which has been used 
in many variants for describing the pseudogap phase of 
the cuprates.^^"^* However, in contrast to these cuprate 
models defined on the square lattice, a gap can be 
opened over the whole Brillouin zone in certain other 
lattices^^"''^'^^ allowing one to characterize the phase by a 
topological invariant. For example, in Ref. [42] a double- 
exchange ferromagnet has been studied on the kagome 
lattice and the ground state has been described as a chi- 
ral spin state with a finite Chern number. Later, Raghu 



et al. studied an extended Hubbard model on the honey- 
comb lattice and showed that both a quantum anomalous 
Hall phase and a quantum spin Hall phase can be gen- 
erated dynamically."^^ A similar idea has also been used 
to obtain a three-dimensional example of an interaction- 
driven topological insulator on the diamond lattice. "^^ 

In our paper, we focus on spinless and spinful ex- 
tended Hubbard models with repulsive interactions on 
the kagome and decorated honeycomb lattice. Inter- 
acting electrons on the kagome lattice provide a model 
system where ferromagnetism can be rigorously shown 
for certain parameters {flat-ban(^^''^'^ and kinetic^^ fer- 
romagnetism). Furthermore, the Mott transition in the 
standard Hubbard model defined on this lattice has been 
studied. ^^ In addition to these examples, a great deal 
of the theoretical work on Hubbard and extended Hub- 
bard models has focused on the case of half-filling where 
the low-energy degrees of freedom are described by a 
frustrated quantum spin model. ^^^^"^ These studies are 
motivated in part by the recent discovery that herbert- 
smithite, a spin-1/2 kagome antiferromagnet, might sup- 
port a spin liquid ground state. ^^'^^ Another system 
where the physics of interacting electrons on the kagome 
lattice might be important is Na2;Co02 where the or- 
bital degrees of freedom give rise to four interpenetrating 
kagome systems.^'^'^'* 

The decorated honeycomb lattice can be viewed as an 
interpolating lattice between honeycomb and kagome. 
While there are few known examples of this lattice in 
nature, ^^ the exact ground states of the Kitaev model 
on this lattice have been found by Yao et.al^^ and other 
higher symmetry spin models have been studied as well.^^ 
Yao et.al^^ has shown that the exact ground state of 
the Kitaev model on this lattice is a chiral spin liq- 
uid that spontaneously breaks time-reversal symmetry. 
There are two topologically distinct chiral spin liquid 
phases: (i) a topologically nontrivial phase with odd 
Chern number and non-abelian vortex excitations and 
(ii) a topologically-trivial phase with even Chern number 
and abelian vortex excitations. In our previous work,^^ 
we have found that this lattice also supports a TI phase 
in the presence of spin-orbit coupling at various filling 
fractions, and we established a connection between the 
topologically nontrivial chiral spin liquid state of the 
Kitaev model (appropriate for strongly interacting elec- 
trons with spin-orbit coupling) and the ground state of 
Z2 topological band insulators (studied in the noninter- 
acting limit). 

Both the kagome and decorated honeycomb lattices 
support a TI in a single-orbital tight-binding model with 
spin-orbit coupling. ^^'^^ In this paper we show that a 
TI (quantum anomalous Hall state for the spinless case) 
can also be interaction-driven on both lattices. We focus 
on filling fractions which either involve a pair of Dirac 
points (1/3 filling in the kagome system) or a quadratic 
band crossing point (2/3 filling in the kagome and 1/2 
filling in the decorated honeycomb system) in the non- 
interacting tight-binding model. ^'^ Using a Hartree-Fock 



mean-field approach, we discuss various possible symme- 
try broken states, present the phase diagrams and high- 
light the competition between various states. We find 
pronounced differences for different filling fractions. In 
particular, a topologically nontrivial phase is the lead- 
ing instability at 2/3 filling on the kagome lattice and 
1/2 filling on the decorated honeycomb lattice. On the 
other hand, to stabilize a topologically nontrivial phase 
at 1/3 filling on the kagome lattice, some fine tuning of 
the interaction parameters is required. We also point out 
that the kagome and decorated honeycomb lattices pro- 
vide examples where topological phases can emerge solely 
due to a complex nearest neighbor hopping, in contrast 
to the honeycomb or diamond lattice in which a complex 
second neighbor hopping is required, at least within a 
single band model. 

The paper is organized as follows. In Sec. II, we intro- 
duce spinless and spinful extended Hubbard models on 
the kagome and decorated honeycomb lattices, and re- 
view the tight-binding band structures and Hartree-Fock 
mean-field approach for the implementation of numeri- 
cal calculations. In Sec. HI and IV, we discuss several 
symmetry-breaking candidate phases and present phase 
diagrams of spinless and spinful extended Hubbard mod- 
els at 1/3 and 2/3 filling fractions. We find the topolog- 
ically nontrivial phases can be stabilized under suitable 
circumstances. Comparisons are also made to related 
work. Then, in Sec. V we briefly discuss the spinless ex- 
tended Hubbard model on the decorated honeycomb lat- 
tice. Finally, we present our conclusions and summary in 
Sec. VI. 



II. MODELS AND METHODS 

We first introduce the models which will be studied 
later by means of the Hartree-Fock approximation. We 
consider both spinless (spin-polarized) and spinful inter- 
acting fermions in a single-orbital Hamiltonian on the 
kagome and the decorated honeycomb lattice. 



A. Extended Hubbard models 

The lattice model under consideration for spinless 
(spin-polarized) fermions takes the form 

(id) (id) 
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Here, c\ annihilates (creates) a spinless fermion on site 
i and Ui = clc^ is the fermion density operator on site 
i. The sums run over nearest- neighbor (i,j), second- 
neighbor ((i,j)), or third-neighbor bonds (((«,j))). The 
hopping amplitude is denoted by t and the parameters 



Vi, V2, and V3 quantify the nearest-neighbor, second- 
neighbor and third-neighbor repulsion, respectively. For 
most parts of our work we set V3 = 0. However, as we 
show later, a small but finite V3 is necessary to stabilize a 
topologically non-trivial insulator for 1/3 filling fraction 
on the kagome lattice. 

The model for spinful fermions includes an additional 
on-site repulsive interaction U. The Hamiltonian reads 

+Vi ^ niUj + ^2^ niUj + V3 ^ niUj. (2) 
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((«J')) 
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„(t) 



Here, cYJ annihilates (creates) a fermion on site i with 

spin a ^t,l, "j<T = cJo-Cio- and Ui = J^cr'^i'^ The sum- 
ming convention and the meaning of the parameters Vi , 
V2, and V3 are the same as for the spinless model. 



B. Kagome and decorated honeycomb lattice 



where fci = ai ■ k and k2 = 0,2 ■ k. There are two dis- 
persing bands (n = 1 and 2) and a flat band (n = 3). At 
fining fraction / — 1/3, the two dispersing bands touch 
at two inequivalent Dirac points located at corners of the 
Brillouin zone 



K± = ±(bi - 62)/3. 



(7) 



At filling fraction / = 2/3, the second band touches the 
flat band at the T point [k = (0, 0)]. This is a quadratic 
band crossing point (QBCP).''^ Upon inclusion of an in- 
trinsic spin-orbit coupling (modeled by a spin-dependent 
imaginary second-neighbor hopping) one finds that a gap 
is opened both at the Dirac points (/ = 1/3) and the 
QBCP (/ = 2/3). 59 The resulting insulating state at 
/ = 1/3 and / = 2/3 is a Z2 topological insulator with 
time-reversal symmetry protected edge states. ^^ In the 
following sections, we explore the possibility of dynam- 
ically generating a topological insulator phase from in- 
teractions and study its competition with other broken- 
symmetry phases. We therefore focus on / = 1/3 and 
/ = 2/3 in this paper. 



The models in Eqs. (1) and (2) have been studied on 
the kagome and the decorated honeycomb lattice in the 
non-interacting limit. ^^'^^ A section of the kagome lattice 
is shown in Fig. 1 and a section of the decorated honey- 
comb lattice is shown in Fig. 11. Both lattices share an 
underlying triangular lattice and we choose the unit cell 
vectors to be 

Qi = (a,0) and 02 = (-, -^a), (3) 

where a is their length. The kagome lattice has three 
sites in the unit cell whereas the decorated honeycomb 
lattice has six. The reciprocal lattice vectors are given 
by 



bi 



27r 



(1, 



and 69 



27r 



(0,^)- 



(4) 



The first Brillouin zone forms a hexagon in momentum 
space for both lattices, similar to the honeycomb lattice 
which also shares the underlying triangular lattice. 



1. Tight-binding band structure on kagome lattice 

The noninteracting energy dispersion for a nearest- 
neighbor tight-binding model [first term in Eq. (1)] can 
be obtained analytically. On the kagome lattice, three 
bands are found with the following dispersion relation: 

ei{k) = -t-tAk, t2{k) = -t + tAk, e3(fc) = 2t. (5) 

In Eq. (5) we have defined 

Af. — \/3 + 2 cos fci + 2 cos ^2 + 2 cos(fci — fc2), (6) 



Tight-binding band structure on decorated honeycomb 
lattice 



Diagonalization of the noninteracting tight-binding 
model on the decorated honeycomb lattice gives the fol- 
lowing six bands: 



ei(fc) = -^- J^t2+t'2+tt'^fc, 



e2{k) = -^--^^-t^+t'^-tt'Au, 
e3{k)=t-t', 

t 19 



eA{k)^-- + s^-t^+t'^-tt'Ak, 

£5(fc) = -\ + ^J^t^+t'^+tt'Ak, 
eeik) =t + t', 



(8a 

(8b 
(8c 

(8d; 
(8e 

(8f: 



and Ak is defined in Eq. (6). Here, we have introduced 
independent hopping amplitudes for hopping within a 
triangle (i) and between triangles (i').^* There are two 
flat bands (n = 3, 6) and four dispersing bands (n = 
1,2,4,5). For filling fractions 1/6 and 2/3, there are 
Dirac points located at K± in the momentum space. 
There are also two quadratic band touching points at 
k = (0,0). The lower QBCP appears at / = 1/2 if 
t' < 3t/2 and at / = 1/3 if t' > 3t/2. The upper QBCP 
appears at / = 5/6. In the presence of a spin-orbit cou- 
pling, TI phases are found at various filling fractions. ^^ 
In this paper we set t = t' and solely focus on / = 1/2. 
Half filling is of particular interest because a topological 
connection between the chiral spin liquid states recently 
found in the Kitaev model'"'^ and the Z2 topological band 
insulator has been established. ^^ 



C. Hartree-Fock mean-field approximation 

We use the standard Hartree-Fock mean-field ap- 
proach to decouple the interaction terms in Eqs. (1) and 
(2). In contrast to comparable studies on the honey- 
comb lattice^^'^^ and the diamond lattice, ^^ we treat the 
Hartree and Fock terms on equal footing in all phases. 



1. Hartree-Fock approximation in the spinless models 

For spinless fermions, we decouple the interaction both 
in the direct and the exchange channel: 



n,n. 



~clcj{c]c,) - {clcj)c]c^ + {4cj){c]c^). (9) 



This procedure yields a mean-field Hamiltonian which is 
bilinear in the fermionic operators and can be diagonal- 
ized. In the following, we focus on uniform phases which 
are characterized by a (possibly enlarged) unit cell. We 
work in the canonical ensemble with a fixed number of 
electrons N^. The free energy at temperature ksT — j3^^ 
is given by 



1 



„/3(Bfc„-p) 



F^-fcsT^log 

.71 



Vi 






+ {V2,V3)-teYms, 



^iN^ 



(10) 
The terms in 



where the chemical potential ^ — fi{T, N^ 
the second and third line of Eq. (10) arise from {ni){nj) 
and {clcj){CjCi) in the decoupling Eq. (9) and are not 
included in the single-particle energies E^n- The most 
general self-consistency (mean-field) equations are 



dF 



OF 



d{ni) d{> 
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OF 

d{c]ci 



= 0. 



(11) 



In the following sections, we discuss various solutions of 
these equations. 



2. Hartree-Fock approximation in the spinful models 

For spinful fermions, we decouple the on-site interac- 
tion according to 

- cltC»;(c4CiT) - {4^C^i)4i(^^^ + (cI^C^) (4iC»t) -(12) 

We assume the mean-field solutions are described by a co- 
linear spin alignment and therefore, without loss of gen- 
erality, we set (cLci^) = {c\\Ci^^) = in what follows. For 
the model on the kagome lattice with Vi = V2 = V3 = 0, 
and at filling fractions / = 1/3 and / = 2/3, we have 
explicitly checked that with all the terms (including 



{<y^^Ci\)) all our self-consistent solutions indeed have a 
co-linear spin alignment. We expect this property will 
persist also for finite further neighbor interactions. How- 
ever, the (cLci^) term has to be kept if one works at 
half filling on the kagome lattice where at the mean-field 
level a coplanar 120° antiferromagnetic state arises in the 
large U limit. The same antiferromagnetic state has also 
been found on the triangular lattice. ^^"^"^ 

The further-neighbor interaction is decoupled in a sim- 
ilar way: 



+ (cLcj73)4^Ci„ - (cLcj/})(c]^c„)j . (13) 

Again, as mentioned above, we set {c[^Cjfj) = for a 7^ /3 
which is justified if the spin alignment is co-linear in the 
physical solutions. The structure of the free energy and 
the self-consistency equations are similar to Eqs. (10) and 
(11) for the spinless models. 



III. SPINLESS FERMIONS ON KAGOME 
LATTICE 

In this section we discuss the zero temperature 
Hartree-Fock mean-field phase diagrams at filling frac- 
tions / ~ 1/3 and / = 2/3 for the spinless model on the 
kagome lattice. We first introduce the candidate phases 
and then show the V1-V2 phase diagrams with and with- 
out a finite V3 for the two special filling fractions. Be- 
cause at / = 1/3 there are Dirac points involved, and at 
/ = 2/3 there is a QBCP, the phase diagrams look rather 
different for these two cases. 



A. Candidate phases 

Let us now introduce possible candidate phases for 
the spinless model. Besides the topologically non-trivial 
quantum anomalous Hall (QAH) phase we also take into 
account possible charge density wave (CDW) patterns. 



1. Quantum anomalous Hall phase 

A complex Fock term in Eq. (9) breaks time-reversal 
symmetry and can give rise to a topological phase charac- 
terized by a non-vanishing Chern number'^^ even though 
there is no external magnetic field. In the present case, 
the total flux through the unit cell must be zero (this 
follows from periodic boundary conditions on the unit 
cell). However, there are finite fluxes through the ele- 
mentary loops and the system shows an integer quantum 
Hall effect. This is in full analogy to Haldane's model on 
the honeycomb lattice.'"' Such a state of matter is called 




QAH / TI 



FIG. 1. (Color online) The spinless flux pattern developed by 
nearest and second nearest interactions that preserves lattice 
symmetry but spontaneously breaks time-reversal symmetry 
on the kagome lattice. Charges are uniform on all sites. The 
blue solid(red dash) line represents a (second)nearest neigh- 
bor complex hopping. For the spinful case, two copies of 
the same(opposite) flux patterns for spin-up and spin-down 
fermions form the quantum anomalous Hall (topological in- 
sulator) state. 




CDWI 



CDWII 



CDwm 



FIG. 2. (Color online) Three different charge density wave 
patterns on the kagome lattice studied in this paper. They 
are characterized by the wave vectors q^ = (0,0), q^ — 62/2 
and qj[[ = (bi — b2)/3. Blue sites stand for the fermion-rich 
(poor) sites at 1/3 (2/3) filling fraction and white sites for 
the fermion-poor (rich) sites at 1/3 (2/3) filling fraction. The 
bond expectation values oscillate in the real space as well. We 
distinguish strong and weak bonds by thick and thin lines. For 
simplicity, we do not show the second neighbor bonds. 



2. Charge density waves 



a quantum anomalous Hall phase and a schematic illus- 
tration of its microscopic current pattern on the kagome 
lattice with finite Vi and V2 is shown in Fig. 1. 

The QAH phase preserves the translational symme- 
try of the noninteracting model but breaks time-reversal 
symmetry. A solution of the self-consistency Eqs. (11) 
is obtained by assuming a uniform charge distribution 
and introducing complex bond expectation values. For 
nearest neighbor bonds we make the following ansatz: 



(4cj) = Xexp(i¥'»j) = Xi + iX2- 



(14) 



A similar ansatz is also made for second-nearest neighbor 
bonds: 



(cicj) = x'exp(i'yJ-,-) = xi + iX2- 



(15) 



There is a gauge freedom in choosing the phase factors 
ipij and Lp[,j because only the inclosed fluxes through el- 
ementary loops are gauge invariant. We choose a uni- 
form gauge 99^' — -inp'^'^ where the sign is fixed accord- 
ing to the dictions of the arrows in Fig. 1. We stress 
that on the kagome lattice a complex nearest-neighbor 
hopping can already stabilize a topologically non-trivial 
phase showing an integer quantum Hall effect. This 
possibility has been explored in a model of a ferromag- 
net with spin anisotropy.*^ Therefore, in contrast to the 
honeycomb'^^'^^ and the diamond lattice, ■^^ the nearest- 
neighbor interaction Vi alone can in principle generate 
a QAH phase if the time reversal symmetry is sponta- 
neously broken. Indeed, we show below that at / = 2/3 
this is the case. However, at f — 1/3 we find it essential 
to have a finite V2 and small V3 in order to stabilize the 
QAH state. 



An effective way to lower the potential energy is to 
develop an inhomogeneous charge distribution. In the 
atomic limit t = and in the absence of further neigh- 
bor interactions, V2 — V3 — 0, there is a macroscopically 
degenerate set of charge configurations which minimize 
the energy. At / = 1/3 (/ = 2/3) these configurations 
obey the "one particle (hole) per triangle"- rule. A finite 
t lifts the degeneracy and in the limit t/Vi ^ 1 the sys- 
tem is effectively described by a hardcore dimer model 
on the honeycomb lattice. ^^'^^'^^ Its ground state is the 
"plaquette" phase with resonating plaquettes and a pe- 
riodicity which triples the unit cell.^^ Physically, it is the 
ring exchange of order |ip/V\^ which stabilizes the pla- 
quette phase. We also note that in the limit t/Vi ^ 1 the 
system becomes particle-hole symmetric. This property 
is clearly lost for small to intermediate interactions, see 
below. 

Further neighbor interactions V2, V3 > lift the degen- 
eracy of the charge configurations in the atomic limit. 
This fact complicates a mapping to an effective dimer 
model for finite t. In the following we study the mean- 
field solutions of a limited number of different classical 
charge distributions. Specifically, we consider three dif- 
ferent charge density wave (CDW) patterns which we de- 
note by I, II and III; see Fig. 2. They were introduced in 
Ref. [65] in order to numerically study the role of the ring 
exchange. For us it is important to realize that V2 > 
favors CDW I as compared to CDW II and III in the 
atomic limit. On the other hand, a third neighbor inter- 
action F3 > favors CDW HI over CDW I and CDW II. 
While the unit cell of pattern I is equal to the nonin- 
teracting unit cell, the unit cell of pattern II is doubled 
and the one of pattern HI is tripled. Note that CDW HI 
can be viewed as the classical charge distribution which 
corresponds to the plaquette phase of the effective dimer 
model in the limit t/Vi ^ 1. 



For the CDW I, the wave vector specifying its peri- 
odicity is Qi = (0, 0) and the densities on the three in- 
equivalent sites of the noninteracting unit ceU are given 
by 



("-l(r„m))/ = f + Pi, 
(n-2(r„m))/ = f + P2, 
(nsirnm))! = f + P3, 



(16) 



where r,„„ = nai + ma2 with (n, m) G Z x Z, / is 
the fiUing fraction and pi + P2 + Ps = 0. Similarly, the 
densities in the CDW II configuration can be written as 



(ni(r'„m))// = / + Pi cos(r, 
{n2{r riTTi)) II = f + P2 cos(r, 
{n3irnm))ii ^ f + Pi cos(r. 
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Qii) 



(17) 



where we have introduced the wave vector qjj — 62/2. In 
our mean-field calculations we find that mirror symmetric 
charge configurations are always favored. Such configura- 
tions are obtained by setting pi = 2p and P2 ^ Ps = —p 
(or cyclically permuted) in Eqs. (16) and (17). CDW I 
and II both break the six- fold rotations symmetry (Cg) 
of the kagome lattice; CDW I breaks it down to C2, while 
CDW II breaks it down even futher. In both cases, there 
are three different possibilities to choose a mirror sym- 
metry plane. The CDW order parameter therefore has 
an additional Z^ freedom. 

Another phase with a mirror symmetric configuration 
of the densities in the non-interacting unit cell is de- 
scribed by the CDW III pattern: 



{ni{rnm))iii = f + 2pcos 

{n2{rnni))lll = / + 2pcos 
{n3{rnm)}lll = / + 2pcos 



27r' 
''nm ■ Qui + s — 

''nm ■ Qui + (s — 1 

"nm ■ Qui + (s + 1 



27r 

T 

2n 

T 
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where the wave vector is Qm = (bi — 62)/3. We have 
introduced the parameter s — 0, 1, 2 which characterizes 
the Z3 freedom in the CDW III. Changing the value of s 
results in a shift of the pattern as a whole either by Oi 
or 02. 

CDW I does not break the original translation sym- 
metry and it can alternatively be viewed as a nematic 
phase. '^^ The direction associated with the nematic order 
is given by 



e = {Q.,Qy)/jQl + Ql 



(19) 



where the components are obtained from the charge and 
bond order 



^x 
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V3 



'^2) _^ o {cW, 



(cU) 
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(20) 



(ni) + (n2)-2(n3) ^^ (clc3) + (4c3)-2(clc2) 



where the "— " sign refers to the case of 1/3 filling 
presently being considered, and the "-I-" sign refers to 
the case of 2/3 filling involving a quadratic band touch- 
ing point. This definition of the nematic order parame- 
ter is in agreement with the definition given in Eq. (4) of 
Ref. [37] for the case of 2/3 filling. 

In our study, we assume that the real hopping expec- 
tation values Xij — {^i^j) obey the same symmetry as 
the charge distribution. In Fig. 2, the weak and strong 
nearest-neighbor bonds are schematically shown. We find 
that taking into account this bond order can significantly 
lower the energy as compared to the case where only the 
Hartree term is kept. For the mirror symmetric solutions, 
the unit vector in Eq. (19) assumes only three different 
directions: 



ei 



(\/3,l)/2, 



62 



(-%/3,l)/2, 63^(0,-1). 



(21) 
These unit vectors will also appear in the low energy 
description of the CDW phases. 



3. Other phases 

Let us now briefly comment on other possible phases 
which are not stabilized in the present models. Dimer- 
ized and trimerized phases were considered in Ref. [59] as 
perturbations to the TI phase in the noninteracting limit. 
We find that for a self-consistency solution with dimer- 
ized or trimerized bonds, it is crucial to take into account 
the charge ordering which results from the bond order. 
However, our numerical results suggest that the charge 
density wave patterns shown in Fig. 2 (where the bond 
order has the same symmetry as the charge order) have 
lower energies than the dimerized or trimerized states. 
We also note that we do not find a mixed QAH and 
CDW phase on kagome lattice, which is in contrast to 
the findings on the checkerboard lattice. ^^ 



B. Phase diagrams at 1/3 filling fraction 

Figure. 3 shows the V1-V2 phase diagrams for (a) 
1/3 = and (b) F3 = 0.4i. At 1/3 fiUing, the noninter- 
acting Fermi "surface" consists of a pair of Dirac points 
located at K± and the density of states vanishes linearly 
at the Fermi energy. As in related studies, '^^'"^^ our mean- 
field calculations yield a stable semi-metallic (SM) phase 
for small to intermediate interactions which can be at- 
tributed to the absence of density of states at the Fermi 
level in the non-interacting limit. 



1. CDW phases and nematic order at f — 1/3 

For large interactions, a CDW phase is stabilized. We 
find CDW I for large Vi and V2 because both the nearest 
neighbor and second neighbor interaction favors CDW I. 
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FIG. 3. The phase diagram of the extended Hubbard model 
for spinless fermions at 1/3 filling fraction on the kagome lat- 
tice. The third neighbor interaction is (a) Vs = and (b) 
V3 = OAt. SM denotes the semi-metallic phase with two 
Dirac points, QAH denotes a time-reversal symmetry broken 
quantum anomalous Hall phase and CDW I and III are charge 
density waves with patterns shown in Fig. 2. Solid lines de- 
note first and dashed lines second order transitions. 



The transition from the SM to the CDW I is first order 
which is different from the situation on the honeycomb 
lattice. -^^ Below we discuss this aspect in more detail. On 
the other hand, CDW III is favored for small V2 and large 
Vi. The transition from the SM to CDW III is second 
order. 

The self-consistent CDW solutions at / = 1/3 are al- 
ways gapped (this is in contrast to f — 2/3 where CDW 
phases with two nodes appear, see Sec. Ill C) and it is 
instructive to look at the corresponding low energy mod- 
els. We first consider the possibility of a weak CDW I 
phase and then argue that it is energetically not favored. 
In fact, only for a large enough order parameter, does the 
CDW I solution have lower energy than the SM phase. 
For simplicity, we keep only the Hartree terms. In lowest 
order in Vp {V — Vi + V2 — 2V3), the effective low energy 
Hamiltonian for the two nodes Z = ± is given by 



Hl=V Y. 4a [^ 



k,l,a,l3 



{k-lA)Uckip+^^\A\\ (22) 



where the velocity is w = \/ita/2 and t = {Tx,Ty) are 
Pauli matrices in the effective "sublattice" space. Fur- 
therrmore, we have introduced an "axial gauge field" ^^ 
A. which can be expressed in terms of the CDW order pa- 
rameter p and the vector e„ specifying the Z3 freedom: 



A= p[nz X e„), 

V 



(23) 



where the e„ are given by Eq. (21). This field shifts the 
position of the Dirac nodes with respect to their orig- 
inal position at K± and consequently, the CDW I de- 
scribed by Eq. (22) has nodes. However, Eq. (22) also in- 
cludes the electron-electron interaction in the mean-field 
description which gives rise to the second term. This 



term is proportional to v'^/V and can be viewed as a 
mass-term for the gauge field. In other words, shifting 
the nodes by the vector A costs an energy proportional 
to l^p. Therefore, it is energetically not favorable to 
built up a finite field A and the SM phase is stable. But 
once V is big enough, the description in terms of Eq. (22) 
breaks down. Solving the full self-consistency equations, 
we find a first order transition from the SM to the gapped 
CDW I phase. 

Let us now consider CDW III which is stable for small 
V2. The wave vector qjjj of pattern III connects the two 
inequivalent Dirac points at K± . From a weak-coupling 
point of view, CDW III therefore opens a gap by cou- 
pling the two Dirac points. This can be made explicit 
by studying the low energy mean-field Bloch Hamilto- 
nian. For simplicity, we set T^2 = V3 = and consider 
only the Hartree terms. The Bloch Hamiltonian for the 
low energy degrees of freedom is expressed in the 4x4 
matrix, 



^™^^)-lv At h(-){k) 
Here, the Dirac Hamiltonians at Kj- are given by 
^(±)(fc) ^vk-T. 



(24) 



(25) 



The coupling between the two Dirac cones can be brought 
into the following form 



2Vip 



1 
-1 



(26) 



Equation (24) can be diagonalized and we find the fol- 
lowing doubly degenerate energy bands 



E±{k)=±Jv^k'^+ApW^ 



(27) 



In particular, an arbitrarily small coupling Vip continu- 
ously opens a gap 41^1 |yo| at the T point in the reduced 
Brillouin zone. (The enlarged unit cell of CDW HI moves 
the low-energy point from Jf-j- to F.) This means that 
the CDW HI is a low-energy instability of the SM phase 
and explains why we observe a second order transition 
at a critical interaction strength. We also note that 
the low-energy theory for the CDW HI, Eq. (24), car- 
ries similarities with the one found for the Kekule tex- 
ture on the honeycomb lattice^^ or in the 7r-flux model 
on the square lattice. ^^ In analogy with these examples, 
we expect that topological defects of the CDW HI pat- 
tern in the form of a Z3 vortex can give rise to interest- 
ing physics; potentially including charge fractionalization 
and anyon statistics. ^^ 



2. Topological phase at f — 1/3 

As shown in Fig. 3(b) we find that a QAH phase can be 
stabilized in a certain region of parameter space. Never- 
theless, it requires some fine tuning of the different inter- 
action strengths. First, we do not find a QAH solution for 
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FIG. 4. (Color online) The phase diagram of the extended 
Hubbard model for spinless fermions at 2/3 filling on the 
kagome lattice. The dotted line indicates where the gap opens 
in the CDW I phase. For V2 = 0, CDW I and II coexist in the 
gray region and CDW III in the black region. We set V3 = 0. 



Vi alone in the parameter space we considered (which is 
different from what we find at 2/3 filling, see Sec. IIIC). 
Second, for a moderate V2 there exists a self-consistent 
solution of Eqs. (11) which breaks time reversal symme- 
try. This QAH phase is triggered by X2 and X21 the 
imaginary part of the nearest and second-neighbor hop- 
ping expectation values (which in general also acquire fi- 
nite real parts). It turns out that for V3 = 0, the CDW I 
phase has lower energy compared to the QAH solution. 
However, a finite V3 increases the energy of the CDW I 
solution making the interaction-driven QAH phase the 
ground state for small Vi and V3 and moderate V2, as 
shown in Fig. 3(b). 



C. Phase diagram at 2/3 filling 

The phase diagram of the extended Hubbard model for 
spinless fermions at 2/3 filling is shown in Fig 4. Here, 
we set V3 = 0. The important difference with 1/3 filling 
is that the Fermi energy in the noninteracting case lies 
at a QBCP between a dispersing and a fiat band. As a 
consequence, the density of states is finite at the Fermi 
energy and the system is unstable to arbitrarily weak 
interactions.'^^ In particular, the semi-metallic phase does 
not survive even for small values of the interactions. The 
phase diagram for low to intermediate interactions looks 
therefore quite different than the corresponding phase 
diagram at 1/3 filling. 



1. Topological phase at f — 2/3 

For small to intermediate interactions we find that the 
QAH phase has the lowest energy. This is in agreement 
with quite general arguments made about the stability 
of a QBCP.'^^ We have numerically calculated the Chern 



number^" associated with this state and found that it is 
±1, indicating it is indeed a topological state displaying 
an integer quantum Hall effect. Note that Vi alone is 
enough to generate the QAH phase because of the par- 
ticular geometry of the kagome lattice with a triangle in 
the unit cell. 

Although the QAH phase is the ground state in a 
rather large region of parameter space, its gap is expo- 
nentially small. The exponential dependence in mean- 
field theory can be found by analyzing the gap equa- 
tion derived from an effective two band Hamiltonian de- 
scribing the low-energy behavior around the QBCP. Let 
us for simplicity set V2 = in the following. A finite 
imaginary part of the nearest-neighbor bond hopping, 
X2 — Iin(CiCj) ^ 0, couples the two bands thereby open- 
ing a gap. In lowest order in Vi, the matrix describing 
this coupling is given by 



ff, 



QAH 



(fc) 



-2iV3I/iX2 



2^^/3ylX2 



(28) 



where £2 zi^) i^ obtained from Eq. (5) by replacing t by 



t + ViXo- 



(29) 



Here, xo = (clcj)o = 1/6 denotes the nearest neighbor 
hopping expectation value in the noninteracting model 
and Eq. (29) takes into account the effect of the Fock 
term in lowest order in Vi . The self-consistency equation 
for X2 reads 



p2t*+K/2 
J2t'-A/2 V(2t* 



N{e) 



2^,2 ' 



isv.'xi 



(30) 



where A is a cutoff energy of the order of t* which is not 
accessible in the low energy description, and iV(e) is the 
noninteracting density of states. Solving Eq. (30) for the 
order parameter X2 yields 



A 



X2 



2^/3yl 



-l/(ViAro) 



(31) 



which holds for small values of the dimensionless coupling 
constant ViNq. Here, we have introduced the density of 
states at the QBCP:^^ 



No = N{2t*) = 



2nt* 



(32) 



The gap is proportional to the order parameter X2 and 
from the result Eq. (31) and the eigenvalues of Eq. (28) 
it follows that 



A 



QAH 



aVsVi 



X2 



2Ae 



-l/(ViN„) 



(33) 



We have checked that the exponential dependence given 
in Eq. (31) is indeed consistent with our full numerical 
evaluation of the self-consistency equations. A similar 
exponential dependence is also found in a one-loop renor- 
malization group treatment'^^ although the dimensionless 
coupling is renormalized compared to Eq. (33). 
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FIG. 5. (Color online) The fluxes $1,2,3 through elementary 
triangles in the QAH phase at fiUing fraction f = 2/3. These 
elementary triangles form the unit cell, as shown in the inset. 
Because of the periodic boundary conditions on the unit cell, 
the net flux is zero and the individual fluxes satisfy 2"l>i + 
$2 + 3$3 = 0. We have set V2 = Vi/2 and V3 = 0. 



We now come back to the general situation where both 
Vi and V2 are finite. In general, the QAH phase is driven 
by both a complex first and second neighbor hopping 
expectation value. Furthermore, one can define an ex- 
plicit deformation^^ of a tight-binding model with com- 
plex nearest neighbor hopping on the kagome lattice to 
show that its ground state is adiabatically connected to 
the ground state of a model with real nearest-neighbor 
hopping and only complex second nearest neighbor hop- 
ping. Therefore, the QAH phase generated by V2 belongs 
to the same topological class as the one generated by Vi . 
Figure 5 shows the fluxes $1,2.3 through three elemen- 
tary triangles forming the unit cell. In this figure, we set 
V2 — 14/2. Because of the periodic boundary conditions 
on a unit ceU, the fluxes satisfy 2$i + <i>2 -f 3$3 = 0. 
Moreover, they are all finite indicating the presence of 
an imaginary hopping amplitude in both the first and 
the second neighbor effective hopping. 



2. CDW phases at 2/3 filling 

For intermediate to large interaction strengths, a CDW 
phase is stable. At V2 = and large 1^, the CDW 111 
phase has the lowest energy. However, the difference in 
energy per site compared to CDW I is only of the or- 
der 10~"^i and becomes smaller the bigger Vi. As a re- 
sult, a very small but finite V2 is sufficient to stabilize 
CDW 1 over CDW III. In contrast to the situation at fill- 
ing fraction / — 1/3, at f — 2/3 CDW 111 can not profit 
from a "nesting" condition. The energy gain compared 
to CDW I is therefore very small. At V2 = 0, a first- 
order phase transition from a QAH state to a CDW 1 




FIG. 6. (Color online) (a) the CDW order parameter p (b) 
the bond order defined as the difference between two nearest 
strong and weak bonds (c) the gap at 2/3 filling fraction on 
the kagome lattice. We have set V2 — t. (d) the sphtting 
of a QBCP(red circle) into two Dirac points(blue cross) for 
gapless CDW 1. 



(11) state takes place at Vi « 1.47i. Numerically, we can 
not resolve any difference in the energy between CDW 1 
and 11 for V2 — 0. Interestingly, there are nodes in the 
CDW I (II) phase where the gap vanishes. The transi- 
tion from the QAH phase to the CDW I state with nodes 
is an example where a transition from a gapped phase 
(QAH) to a gapless phase (CDW I) occurs by increasing 
the interaction strength, [see Fig. 6(c)]. 

The gapless CDW I persists even for finite V2. Figure 6 
shows the CDW order parameter p, the bond order i' = 
Xs ~ Xw defined as the difference between two nearest 
strong and weak bonds and the magnitude of the gap 
as a function of Vi for fixed V2 = t. We can see that a 
finite CDW order is accompanied by a finite bond order. 
In fact, both types of orders jump to a finite value at 
the transition Vi^ from the QAH to the CDW I phase. 
Note, however, that the gap is zero up to a second critical 
interaction strength Vic2 indicating the presence of band 
degeneracy points below Vic2 ■ At Vic2 a kink is observed 
in the order parameters and the gap gradually starts to 
increase. 

In the following, we show that the gapless CDW 1 phase 
results from the splitting of the QBCP into two nodes. '^^ 
We notice that the change of bond-order is one order of 
magnitude smaller than t and therefore can be neglected 
for the moment. We find the following low-energy Bloch 
Hamiltonian 



H, 



CDW 



(fc) = 



- "=2 



l{k)+t*ucos{2(l)) 



K^ sin(20) 
^sin(20) el{k)-t*ucos{2(l>)^ 

(34) 
where u — 2(14 + V2){p + iv)/t* . The renormalized 
hopping t* is given by Eq. (29) with xo replaced by 
X = {Xs + 2xid)/3 and we have neglected the Fock-terms 
generated by the V2 interaction (this term is negligible in 
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practice). In Eq. (34) we have introduced polar coordi- 
nates {k,(j)) which are defined by fc • e„ = kcoscj). Note 
that right at the F-point the angle cf> is not well-defined. 
Equation (34) should be contrasted with Eq. (28) for the 
QAH phase: as opposed to the QAH order parameter, 
the CDW order parameter introduces an anisotropic an- 
gle dependent effective coupling between the two bands. 
Expanding the dispersion around the T point we find for 
the eigenvalues of Eq. (34), 



E2{k, 4>)/t* = 2 - [fc2 + ^B{k, 0)]/8, 


(35) 


Ez{k, cp)/t* = 2 - [fc2 _ ^B{k, 0)]/8. 


(36) 


function -B(fc, (j)) is given by 




B(fc,0) = k'^ - l&k^ucos{2(j}) + Mu^ . 


(37) 



It has roots at two points where the two bands touches: 



ku = V8u, 



= 0,7r. 



(38) 



This analysis shows that a finite CDW I order splits the 
QBCP into two nodes moving along the line defined by 
the vector e„. The bottom right panel of Fig. 6 illustrates 
the situation for 63. We have calculated the Berry phase 
(winding number)'^^ of the QBCP and found that it is 27r 
{— mod 27r). The corresponding Berry phases (wind- 
ing numbers) of the two nodes appearing in the gapless 
CDW I are both tt. Thus, the QBCP does sphts into 
two Dirac points with Berry phases tt conserving to total 
winding number, as it was suggested in Ref. [37]. 



D. Comparison with existing work 

Recently, several numerical works^^'^^ appeared deal- 
ing with the charge density wave order on the kagome 
lattice dX J — 1/3 or / = 2/3. Here, we want to 
briefly relate our results with their findings. In Nishi- 
moto ei.aL's work,^^ the authors considered the large Vi 
limit with vanishing V2. They showed that CDW III is 
the ground state that is consistent with the "plaquette" 
state obtained from an effective quantum dimer model 
on the honeycomb lattice.^' In this strong interacting 
limit, / = 1/3 and / = 2/3 are equivalent and numeri- 
cal calculations^^ confirm that CDW III is stabilized by 
the ring exchange process proportional to |tp/y^. Inter- 
estingly, fractionalized excitations with charge e/2 have 
recently also been reported in the strong coupling limit. ^^ 

Our mean-field calculation cannot capture the resonat- 
ing nature of the quantum dimer model and is not valid 
in the strongly interacting limit. However, CDW III in 
the mean-field treatment can be viewed as the "classical" 
configuration of plaquette states. At 1/3 filling fraction, 
the CDW III is found to be more stable than either CDW 
II or I at large Fi, and the energy difference between 
them becomes smaller as Vi grows. This is consistent 
with Nishimoto et.al.^s work. Furthermore, we predict a 
metal- insulator transition takes places at Vic = 3. It at 



1/3 filling, which is in quite good agreement with their 
result Vic = 4.0t. However, at 2/3 filling, our mean- 
field results differ significantly from theirs. We find that 
the leading instability at small interactions is the QAH 
state that spontaneously breaks time-reversal symmetry 
and has an exponentially small gap.'^^ A metal-insulator 
transition takes place around Vic = 2t in our study while 
Nishimoto et.al. reported a metal-insulator transition at 
finite Vic = 2.6t. 



IV. SPINFUL MODEL ON THE KAGOME 
LATTICE 

Let us now turn to the spinful model on the kagome 
lattice. The additional spin degrees of freedom add con- 
siderable complexity to the problem and introduce sev- 
eral more potential phases. 



A. Candidate phases 

The SM and the CDW phase are equivalent to those 
in the spinless model. Here, we discuss additional phases 
which appear in the spinful model. 



1. Topological insulator and quantum anomalous Hall state 

The topological insulator and quantum anomalous Hall 
state are both stabilized by a complex Fock term of near- 
est or second neighbor interaction which gives rise to a 
complex hopping amplitude. The difference between TI 
and QAH phases can be described by the 2x2 matrix 
(Ci^Cjp) defined in the spin space as discussed below. 

The QAH state breaks the time reversal symmetry but 
not the spin rotation SU(2) symmetry. Therefore, the 
most general form of the uniform phase consistent with 
these requirements is 



(cLS/s) = [(Xl+«X2)cro]„^ 



(39) 



where Xi ^-^d X2 are real numbers and CTq = 1 is the iden- 
tity matrix. A phase with a finite X2 shows an anomalous 
quantum Hall effect and a non-zero Chern number. 

On the other hand, the TI does not break time reversal 
symmetry but breaks the SU(2) spin rotation symmetry 
down to U(l). The most general form is therefore 



(cLS/j) = ki'^o + *^2(rt • ff)]cp 



(40) 



where xi and X2 are both real numbers and n is a unit 
vector describing how the SU(2) spin rotation symme- 
try is broken. In other words, spin-rotation symmetry 
is only preserved for rotations around n. Without loss 
of generality we can assume n = fiz- We note that al- 
lowing xii X2 and n to be spatially dependent allows 
one to study topological defects of the order parameter. 
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such as skyrmions, providing a potential route for ex- 
otic superconductivity/^ On the other hand, in contrast 
to their two-dimensional counterpart, three-dimensional 
interaction-driven Tls completely break the spin-rotation 
symmetry and their order parameter involves a rotation 
matrix. Again, it is possible to study topological defects 
which host protected modes. ^^ 

A short inspection of the mean-field free energy of 
the TI and the QAH phase shows that these two phases 
are degenerate on the mean-field level. It is likely that 
fiuctuations around the mean-field state might favor one 
phase over the other. Because the TI breaks the contin- 
uous spin-rotation symmetry, there are Goldstone modes 
in the ordered phase. "^^ It was suggested^^ that quan- 
tum fluctuations associated with these modes lower the 
ground state energy of the TI as compared with the QAH 
phase which does not have Goldstone modes. This argu- 
ment appears to be confirmed via "unbiased" functional 
renormalization group methods. '^^ We don't see any rea- 
son for those arguments not to hold in the present case 
as well. 
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FIG. 7. (Color online) Schematic of four types of candi- 
date phases on the kagome lattice for 1/3 filling fraction; 
(a) SCDW I (b) SCDW II (c) bond-order wave (BOW) and 
(d) ferromagnet (FM). Upward arrows and downward arrows 
denote the magnetization on each site. The same/different 
circles represent same/different numbers of fermions on cor- 
responding sites. For simplicity we only show nearest bonds 
(the addition of two spin-resolved bonds) and do not show 
the second nearest bonds. Stronger bonds are shown in bold. 



There is another class of phases that emerges as a result 
of the special filling fractions, the non-bipartite nature 
of the kagome lattice and the additional spin degrees of 
freedom. We term it "spin charge density wave" (SCDW) 
because it involves both a spin and a charge density wave. 
In our mean-field calculations we restrict to phases which 
do not break the translational symmetry. By solving self- 
consistency equations we identify two types of SCDWs 
which are stable for some interaction parameters. 

The first pattern, SCDW I, is characterized by the fol- 
lowing distribution: 



the following distribution: 



("n) = 


= f + p + m, 


(41a) 


("u) = 


= f + p-m, 


(41b) 


("2t> = 


= f + p-m, 


(41c) 


("2i> = 


= f + p + m, 


(41d) 


("3T> = 


= I-^P, 


(41e) 


("3;) = 


= f-'2p. 


(41f) 



Here, p and m are the charge density and spin density 
order parameter, respectively. Furthermore, we assume 
that the symmetry of the bond expectation values is de- 
termined by the symmetry of the spin-charge configu- 
ration and therefore, three different spin-resolved bond 
expectation values have to be introduced. The phase 
SCDW I is schematically shown in Fig. 7(a). 

The other configuration, SCDW II, is characterized by 



("It) = f - P + m, 
("u) = f - p-m, 
("-2t> ^ f - P + m, 
{n2i) ^ f - p-m, 
("3t> = / + 2p - 2m, 
("3;) = / + 2p + 2m. 



(42a) 
(42b) 
(42c) 
(42d) 
(42e) 
(42f) 



The schematics of the SCDW II is shown in Fig. 7(b). In 
addition, we also introduce four different spin-resolved 
bond expectation values to make it consistent with the 
above spin-charge distribution. 

Both SCDWs have zero magnetization in the unit cell. 
However, they differ in that SCDW I has antiferromag- 
netic order in the a^ direction, ferromagnetic order in the 
a2 and a^ — a-2, directions, while SCDW II has ferromag- 
netic ordering in the ax direction but antiferromagnetic 
ordering in the a-i and a^ — 02 directions. 

In the next section, we will see that SCDWs arise in the 
case of large U but small or moderate V\ and V2. At 1/3 
filling fraction, SCDWs can be understood as a means to 
reduce the on-site interactions by single occupancy at two 
sites in a unit cell. Therefore they become unstable when 
V\ or Vi becomes large and CDW dominates. We stress 
that the solutions of SCDWs are saddle points of the free 
energy instead of a global minimum in the usual situa- 
tion, therefore, one has to solve self-consistency equations 
directly to obtain the SCDW solutions. 



12 



3. Bond-order wave 

Next we consider the bond-order wave (BOW) as has 
been found in Ref. [72] for the t-J model at / — 1/3 
under quite general conditions. The BOW is character- 
ized by a uniform charge distribution and a bond or- 
der which breaks the inversion symmetry of the unit cell 
by establishing strong bonds Xs for the up triangles and 
weak bonds for the down triangles Xw ■ It is schematically 
shown in Fig. 7(c). 



4- Ferromagnet 

The ferromagnetic state (FM) is characterized by a 
uniform magnetization density m. The spin densities are 
given by 



{n-it) = f + 171/2, {n,i) = / - m/2, 



and we introduce the Fock terms 



Xt = (cltSt)' Xi = (^US;) 



(43) 



(44) 



The Fock terms are different for nearest-neighbor and 
second-neighbor bonds. Equations (43) and (44) are used 
as an ansatz to solve the self-consistency equations nu- 
merically. A finite magnetization of the form Eq. (43) 
introduces a Zeeman field which uniformly lowers the 
energy of the spin-f electrons with respect to the spin- J, 
electrons by Um. At / = 1/3 the maximally polarized 
state is obtained when there are two up electrons per 
unit cell. At / = 2/3, the maximally polarized state 
corresponds to 3 up electrons and one down electron per 
unit cell. At both filling fractions the saturated value of 
the magnetization is TTigat = 2/3. In the next section, we 
will see that the maximally polarized FM state arises in 
the large U limit. 



B. Phase diagrams at 1/3 filling 

We first discuss the role of U and Vi and set V2 = 
V3 = 0. This allows for a direct comparison with the 
phase diagram at 2/3 filling shown below. 



1. U-Vi phase diagram at 1/3 filling 

The C/-Vi -phase diagram is shown in Fig. 8. Similar to 
the spinless model, we find that the SM is stable for small 
to intermediate interactions which we again attribute to 
the vanishing density of states at the Fermi energy in the 
noninteracting limit. For dominant Vi interaction, we 
find that the CDW III is stable and the transition from 
SM to the CDW III is second order. For dominant onsite 
interaction U, a SCDW phase is stabilized. Both patterns 
SCDW I and SCDW II are stable for some values of the 



b 



12 
10 

8 

6: 

4 
2 







SCDW I 




^^ 


. SCDW II 






30W 


-__/: 


--- " 




J^I 




SM 


,.'' 




CDW III 






.-'' 




■ 



3 

Vi/t 



FIG. 8. The phase diagram of the spinful model at 1/3 filling 
on the kagome lattice. The SCDW 1 and II phase involve 
both a finite charge and spin density wave order parameter. 
Furthermore, when U competes with Vi a bond-order wave 
(BOW) is found. Solid lines indicate first order and dashed 
lines second order transitions. 



interaction. We note that for small Vi there is a second 
order transition from the SM to the SCDW II with nodes. 

For intermediate U and Vi we find BOW is the fa- 
vored ground state on kagome lattice at 1/3 filling. ^^ It 
requires that U is of the same order as Vi to suppress 
the CDW III. On the other hand, it requires a reasonable 
value of Vi to generate the bond order at all. However, we 
expect that the superexchange mechanism (second order 
in t/U), which is not captured in our mean-field treat- 
ment, could stabilize this phase also for smaller ViJ'^ 

At quite large onsite interactions {U ~ 20i) a FM 
phase is stabilized (not shown). The FM state is fully 
polarized at the mean-field level and has an energy gain 
of 



epM — esM — £2 — Ci 



V 

24i 



(e1-e1-2e-ie-2)-| (45) 



per unit cell as compared to the SM phase. In Eq. (45), 
we have introduced the average kinetic energy of the filled 
band n, 



k 



.(fc), 



(46) 



where N is the number of unit cells in the lattice, and 
the dispersion relation e„(fc) is given in Eq. (5). We note 
that the presence of a FM state for large interactions 
is consistent with numerical studies.'*^ Finally, we note 
that the QAH/TI phase does not occur in the absent of 
a finite V2. Again, this is in agreement with the spinless 
case. 



2. U-V2 phase diagram at 1/3 filling 

The U-V2 phase diagram of is shown in Fig. 9. Like in 
the spinless case at 1/3 filling fraction, we add a small 
V3 interaction to suppress CDW I for finite V2 and stabi- 
lize TI/QAH. The overall structure is quite similar to the 
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FIG. 9. The U-V2 phase diagram for Vi = and Vs = 0.4i. An 
interaction-driven TI appears for finite U and V2. Sohd lines 
indicate first order and dashed fines second order transitions 




FIG. 10. The U - Vi phase diagram for ^2 = V3 = at 
2/3 fiUing fraction. Similar to the spinless case at 2/3 filling 
fraction, CDW I has nodes which separates itself from gapped 
phase by a dash-dot line. Similarly, SCDW I has nodes for 
small interaction strengths. 



U-Vi phase diagram. However, the charge density wave 
has pattern I for large V2 since large V2 does not favor 
CDW III but CDW I. BOW phase is now replaced by the 
QAH/TI phase. That the topological phase appears in 
the middle of the phase space seems to be a rather univer- 
sal feature in systems which have a Dirac point and has 
also been reported on the honeycomb and the diamond 
lattice.35.36 poj. j^rge U ~ 6t, SCDW II is stabilized and 
we find that it is gapless. A first-order phase transition 
from SCDW II to SCDW I occurs when V2 increases and 
finally CDW I dominates for large Vi. Note a FM state 
occurs at even larger U ~ 20i (not shown) . 



C. Phase diagram at 2/3 filling 

For 2/3 filling we focus only on the U-Vi phase- 
diagram. The phase diagram is shown in Fig. 10. Most 
importantly, we found that the dominant instability for 
arbitrarily small Vi is to the QAH/TI phase and this 



phase survives also for finite U. Increasing Vi further, 
there is a first order transition to a gapless and then 
gapped CDW. 

We note that for Vi < 0.3t the energy difference be- 
tween various states is very small: SCDW I and II as well 
as QAH/TI have energy differences of less than 10~®i per 
unit cell and we had to use a very high precision in the 
numerical calculation to resolve the phase diagram. How- 
ever, for larger values of U the FM phase is clearly favored 
in the mean-field calculation. This is again a maximally 
polarized FM state. The energy density as compared to 
the SM phase is 
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U 
J' 



(47) 



and we have used the definition Eq. (46). While the 
energy gain for the onsite repulsion U is the same in the 
SCDWs, it is the kinetic energy which favors the FM 
phase over the SCDWs for large U. 



V. SPINLESS FERMIONS ON DECORATED 
HONEYCOMB LATTICE 

In this section, we briefly examine the possibility of 
an interaction-driven QAH state for spinless fermions on 
the decorated honeycomb lattice and discuss the rela- 
tionship of the QAH phase with other competing phases. 
The study of interaction effects on this lattice is par- 
tially motivated by a recent paper that exactly solved^^ 
the Kitaev model on this lattice in the strongly interact- 
ing limit of the underlying fermions; our previous paper 
established the existence of TIs on this lattice in the non- 
interacting limit. ^'^ One natural question to ask is what 
will happen for intermediate interaction strengths where 
the Hartree-Fock mean-field approximation is still valid. 
We work at half filling and t' = t and show that the QAH 
state is the leading instability in the presence of interac- 
tions. Moreover, it occupies a rather wide region in the 
phase diagram. 

The Fermi surface at 1/2 filling lies at a quadratic band 
crossing point in the center of the Brillouin zone, where 
a flat band crosses a quadratic band. It allows the emer- 
gence of a QAH phase quite easily without any fine tun- 
ing of interaction strengths. We consider the nearest- 
neighbor interaction Vi and second nearest-neighbor in- 
teraction V2 on this lattice. For the Vi interaction, we can 
introduce a dynamically generated flux pattern in the two 
triangles. We also introduce a second-neighbor flux in the 
same way as in the kagome lattice [see Fig. 11(a)]. One 
key difference, however, is that we have to allow the possi- 
bility of different values of inter-triangle complex hopping 
parameters and intra-triangle complex hopping parame- 
ters due to nonequivalence of the two hopping parameters 
in the non- interacting limit. One can easily show that if 
the flux through a unit cell is zero and time-reversal sym- 
metry is broken, a QAH state is realized similar to the 
one on kagome lattice. In our calculation, we also find it 
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FIG. 11. (Color online) (a) the flux pattern developed by 
interactions that preserves the lattice symmetry but sponta- 
neously breaks time reversal symmetry on decorated honey- 
comb ( "star" ) lattice. The blue solid line with an arrow (two 
arrows) represents a nearest neighbor intra-triangle(inter- 
triangle) complex hopping while the red dash line with an 
arrow represents the second neighbor complex hopping, (b) 
the favorable CDW pattern from solutions of self-consistency 
equations. We have used different combinations of color and 
markers to show the mirror symmetry. Real first neighbor 
bonds consistent with symmetry of the CDW pattern have 
been assumed (second neighbor bonds not shown). 



is possible to have a BOW state if the phase of the flux is 
zero or tt. The BOW is very close in energy to the QAH 
state (lO^^i), but appears to lose out for the parameter 
ranges we studied. 

We vifill restrict ourselves to q = CDW that origi- 
nates from the Hartree term of the mean-field Hamilto- 
nian. At 1/2 filling, one can see that Vi and V2 frustrate 
each other, in contrast to the kagome lattice at 1/3 filling 
where Vi and V2 both stabifize a CDW state. The CDW 
pattern we found based on the mean-field self-consistency 
equations is shown in Fig 11(b). Among several CDW so- 
lutions we have found, we identify this particular CDW 
with a mirror symmetry as having the lowest ground 
state energy. Real bond orders from the Fock term have 
also been introduced implicitly to be consistent with this 
CDW pattern. 

The phase diagram is shown in Fig. 12. Along the 
horizontal axis where V2 = 0, a BOW phase competes 
with the QAH for small interaction but is higher in en- 
ergy by ~ IQ-^t. A CDW phase develops for Vi > 1.9i. 
Generally, the energy gain from forming a uniform distri- 
bution of fermions and complex bonds between neighbor 
sites is very small compared to that in the CDW phases. 
Therefore, the QAH phase is the only favorable ground 
state when either (i) there are no CDW solutions, or (ii) 
the CDW states are frustrated or suppressed and there- 
fore have much higher ground state energy. The first 
case has been seen in the kagome lattice at 2/3 filling 
fraction for small interactions, and the second situation 
is realized in the present case of the decorated honey- 
comb lattice at 1/2 filling where a QAH phase occupies 
a large region of the phase diagram. This is intimately 




FIG. 12. (Color online) The phase diagram of the extended 
Hubbard model of spinless fermions at 1/2 flUing on the dec- 
orated honeycomb lattice. Due to the mutual frustration of 
Vi and V2 , the QAH phase occupies the middle part of phase 
diagram and the regions of CDW phase have been split into 
two parts. We have set t' = t. 



related to the mutual frustration between the Vi and V2 
interaction at 1/2 filling. To see this, let us consider the 
large V2 S> Vi limit. One finds that in order to mini- 
mize the interaction V2, the preferred ground state is the 
CDW configuration where three sites in a top triangle 
are almost occupied while the three sites in the bottom 
triangle in the same unit cell are almost empty. This is 
exactly the sublattice potential perturbation considered 
in Ref. [58] that destabilizes the quantum spin Hall phase. 
Though this sublattice potential appears artificial at first 
sight, we show here that it can result from a many-body 
interaction. Clearly this configuration is not stable if Vi 
is increased beyond a critical value. This explains the 
fact that at large Vi or V2 the CDW phase is the ground 
state, while the QAH state is the ground state when Vi 
is comparable to V2. 

It is possible to perform a similar mean-field calcula- 
tion for the spinful case, and one expects that a TI/QAH 
phase will dominate at small interactions strengths for 
1/2 filling. However, the details are beyond the scope of 
this paper and is left to future work. 



VI. CONCLUSIONS AND SUMMARY 

We have presented comprehensive Hartree-Fock mean- 
field calculations of the phase diagram for spinless and 
spinful fermions described by the extended Hubbard 
model on the kagome lattice and decorated honeycomb 
lattice. We have established the existence of interaction- 
driven topological phases at filling fractions where either 
Dirac points or quadratic band crossing points are in- 
volved. We find that both TI and QAH phases can be 
described by conventional complex bond order parame- 
ters. Quite generally, we find that at 2/3 filling on the 
kagome lattice and 1/2 filing on the decorated honey- 
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comb lattice (where a quadratic band crossing point is 
involved in the non-interacting limit), the TI/QAH phase 
is the leading instability for small interaction strengths. 
We have observed also that interaction-driven topologi- 
cal phases only exist beyond a critical interaction value 
when the Fermi surface lies at Dirac points at 1/3 filling 
on kagome lattice (in the zero interaction limit). Fur- 
thermore, we discuss in detail various other phases in- 
cluding charge density wave, spin-charge density wave, 
bond-ordered wave, and ferromagnets on the two lattices. 
An important lesson drawn from this study is that 
systems whose non-interacting band structures involve 
quadratic band crossing points can be unstable to topo- 
logical phases with arbitrarily weak interactions, even in 
the absence of microscopic spin-orbit coupling. We hope 



this work will aid in the search for topological states 
of matter by enlarging the class of candidate materials 
to include those which do not have strong microscopic 
spin-orbit coupling but do have certain features (such as 
quadratic band touching points) in their non-interacting 
band structure. 



ACKNOWLEDGMENTS 

We appreciate stimulating discussions with Marcel 
Franz, Mike Hermele, Srinivas Raghu, Kai Sun, and Fan 
Zhang. We gratefully acknowledge financial support from 
ARO Grant W911NF-09-1-0527. 



^ X.-G. Wen, Quantum Field Theory of Many-Body Systems 

(Oxford, New York, 2004). 
^ C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das 

Sarma, Rev. Mod. Phys. 80, 1083 (2008). 
^ G. E. Volovik, Ttie Universe in a Helium Droplet (Claren- 
don Press, Oxford, 2003). 
" D. J. Thouless, M. Kohmoto, M. P. Nightingale, and 

M. den Nijs, Phys. Rev. Lett. 49, 405 (Aug 1982). 
^ F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (Oct 1988). 
^ A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Lud- 

wig, Phys. Rev. B 78, 195125 (Nov 2008). 
^ J. E. Moore, Nature 464, 194 (2010). 
® C. L. Kane and M. Hasan, arXiv:1002.3895vl(2010). 
^ X.-L. Qi and S.-C. Zhang, Physics Today 63, 33 (2010). 
^° C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 

(Nov 2005). 
" C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 

(Sep 2005). 
^^ L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 

106803 (Mar 2007). 
^^ J. E. Moore and L. Balents, Phys. Rev. B 75, 121306 (Mar 

2007). 
^'^ R. Roy Pliys. Rev. B 79, 195322 (May 2009). 
^^ B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Science 

314, 1757 (2006). 
^^ B. A. Bernevig and S.-C. Zhang, Phys. Rev. Lett. 96, 

106802 (2006). 
^'' S. Murakami, S. Iso, Y. Avishai, M. Onoda, and N. Na- 

gaosa, Phys. Rev. B 76, 205304 (2007). 
^* M. Konig, S. Wiedmann, C. Brune, A. Roth, H. Buhmann, 
L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, Science 318, 
766 (2007). 
^^ A. Roth, C. Brune, H. Buhmann, L. W. Molenkamp, 
J. Maciejko, X.-L. Qi, and S.-C. Zhang, Science 325, 294 
(2009). 
2° D. Hsieh, D. Qian, L. Wray Y. Xia, Y. S. Hor, R. J. Cava, 

and M. Z. Hasan, Nature 452, 970 (2008). 
2^ D. Hsieh, Y. Xia, L. Wray D. Qian, A. Pal, J. H. Dil, 
J. Osterwalder, F. Meier, G. Bihlmayer, C. L. Kane, Y. S. 
Hor, R. J. Cava, and M. Z. Hasan, Science 323, 919 (2009). 
^^ Y. Xia, D. Qian, D. Hsieh, L. Wray A. Pal, H. Lin, A. Ban- 
sil, D. Grauer, Y. S. Hor, R. J. Cava, and M. Z. Hasan, 
Nat. Phys. 5, 398 (2009). 



2^ Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S.-K. Mo, 
X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, S. C. 
Zhang, L R. Fisher, Z. Hussain, and Z.-X. Shen, Science 
325, 178 (2009). 

2^ D. Hsieh, Y. Xia, D. Qian, L. Wray J. H. Dil, F. Meier, 
J. Osterwalder, L. Patthey, J. G. Checkelsky, N. P. Ong, 
A. V. Fedorov, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, 
R. J. Cava, and M. Z. Hasan, Nature 460, 1101 (2009). 

25 D. Hsieh, Y. Xia, D. Qian, L. Wray F. Meier, J. H. Dil, 
J. Osterwalder, L. Patthey, A. V. Fedorov, H. Lin, A. Ban- 
sil, D. Grauer, Y. S. Hor, R. J. Cava, and M. Z. Hasan, 
Phys. Rev. Lett. 103, 146401 (2009). 

2« H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. 
Zhang, Nat. Phys. 5, 438 (2009). 

^■^ B. Yan, C.-X. Liu, H.-J. Zhang, X.-L. Q. Chi- Yung Yam, 
T. Frauenheim, and S.-C. Zhang, arXiv: 1003.0074(2010). 

2* S. Chadov, X.-L. Qi, J. Kiibler, G. H. Fecher, C. Felser, 
and S. C. Zhang, arXiv: 1003.0193(2010). 

2® H. Lin, L. Wray, Y. Xia, S. Jia, R. Cava, A. Bansil, and 
M. Hasan, arXiv:1003.0155(2010). 

^° H. Lin, L. Wray Y. Xia, S.-Y. Xu, S. Jia, R. Cava, A. Ban- 
sil, and M. Hasan, arXiv:1004.0999(2010). 

^^ J. Li, R.-L. Chu, J. K. Jain, and S.-Q. Shen, Phys. Rev. 
Lett. 102, 136806 (Apr 2009). 

^2 C. W. Groth, M. Wimmer, A. R. Akhmerov, J. Tworzydlo, 
and C. W. J. Beenakker, Phys. Rev. Lett. 103, 196805 
(Nov 2009). 

^^ D. Pesin and L. Balents, Nat. Phys. 6, 376 (2010). 

3'' S. Rachel and K. L. Hur, arXiv:1003.2238v3(2010). 

^5 S. Raghu, X.-L. Qi, C. Honerkamp, and S.-C. Zhang, Phys. 
Rev. Lett. 100, 156401 (2008). 

^^ Y. Zhang, Y. Ran, and A. Vishwanath, Phys. Rev. B 79, 
245331 (Jun 2009). 

^■^ K. Sun, H. Yao, E. Fradkin, and S. A. Kivelson, Phys. Rev. 
Lett. 103, 046811 (2009). 

^* L Affleck and J. B. Marston, Phys. Rev. B 37, 3774 (Mar 
1988). 

^'^ X.-G. Wen and P. A. Lee, Phys. Rev. Lett. 76, 503 (1996). 

^° C. M. Varma, Phys. Rev. Lett. 83, 3538 (1999). 

^^ S. Chakravarty, R. B. Laughlin, D. K. Morr, and C. Nayak, 
Phys. Rev. B 63, 094503 (2001). 

'^^ K. Ohgushi, S. Murakami, and N. Nagaosa, Phys. Rev. B 
62, R6065 (Sep 2000). 



16 



*^ A. Mielke, J. Pliys. A 25, 4335 (1992). 

'*'' H. Tasaki, Progress of Theoretical Physics 99, 489 (1998). 

*^ F. PoUmann, P. Fulde, and K. Shtengel, Phys. Rev. Lett. 

100, 136404 (Apr 2008). 
'^^ T. Ohashi, N. Kawakami, and H. Tsunetsugu, Phys. Rev. 

Lett. 97, 066401 (Aug 2006). 
^'^ M. Hermele, Y. Ran, P. A. Lee, and X.-G. Wen, Phys. Rev. 

B 77, 224413 (Jun 2008). 
■** F. Wang and A. Vishwanath, Phys. Rev. B 74, 174423 

(2006). 
^^ R. R. P. Singh and D. A. Huse, Phys. Rev. B 76, 180407 

(2007). 
^° G. Evenbly and G. Vidal, Phys. Rev. Lett. 104, 187203 

(2010). 
5^ J. S. Helton, K. Matan, M. P. Shores, E. A. Nytko, B. M. 

Bartlett, Y. Yoshida, Y. Takano, A. Suslov, Y. Qiu, J.-H. 

Chung, D. G. Nocera, and Y. S. Lee, Phys. Rev. Lett. 98, 

107204 (2007). 
^^ T. Imai, E. A. Nytko, B. M. Bartlett, M. P. Shores, and 

D. G. Nocera, Phys. Rev. Lett. 100, 077203 (Feb 2008). 
^^ W. Koshibae and S. Maekawa, Phys. Rev. Lett. 91, 257003 

(Dec 2003). 
^^ M. Indergand, Y. Yamashita, H. Kusunose, and M. Sigrist, 

Phys. Rev. B 71, 214414 (Jun 2005). 
^^ Y.-Z. Zheng, M.-L. Tong, W. Xue, W.-Z. Zhang, X.-M. 

Ghen, F. Grandjean, and G. Long, Agnew. Chem. Int. Edn. 

46, 6076 (2007). 
^^ H. Yao and S. A. Kivelson, Phys. Rev. Lett. 99, 247203 

(2007). 



^'^ B.-J. Yang, A. Paramekanti, and Y. B. Kim, Phys. Rev. 

81, 134418 (2010). 
^^ A. Riiegg, J. Wen, and G. A. Fiete, Phys. Rev. B 81, 

205115 (2010). 
^^ H.-M. Guo and M. Franz, Phys. Rev. B 80, 113102 (2009). 
•^^ M. Kargarian and G. A. Fiete, arXiv:1005.44357(2010). 
•^^ G. Weeks and M. Franz, Phys. Rev. B 81, 085105 (Feb 

2010). 
^^ D. A. Huse and V. Elser, Phys. Rev. Lett, 60, 2531 (Jun 

1988). 
®'^ B. Bernu, P. Lecheminant, G. Lhuillier, and L. Pierre, 

Phys. Rev. B 50, 10048 (Oct 1994). 
•^^ A. Singh, Phys. Rev. B 71, 214406 (Jun 2005). 
'^^ S. Nishimoto, M. Nakamura, A. O'Brien, and P. Fulde, 

arXiv: 1001. 0808(2010). 
^"^ A. O'Brien, F. PoUmann, and P. Fulde, 

arXiv:1002.4894(2010). 
^■^ R. Moessner, S. L. Sondhi, and P. Ghandra, Phys. Rev. B 

64, 144416 (Sep 2001). 
^^ G.-Y. Hon, G. Ghamon, and G. Mudry, Phys. Rev. Lett. 

98, 186809 (May 2007). 
^^ B. Seradjeh and M. Franz, Phys. Rev. Lett. 101, 146401 

(Sep 2008). 
™ T. Fukui, Y. Hatsugai, and H. Suzuki, J. Phys. Soc. Jpn. 

74, 1674 (2005). 
^^ T. Grover and T. Senthil, Phys. Rev. Lett. 100, 156804 

(2008). 
'^^ M. Indergand, A. Lauchli, S. Capponi, and M. Sigrist, 

Phys. Rev. B 74, 064429 (Aug 2006). 



